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ABSTRACT 

Aims. To predict time delays for a sample of gravitationally lensed quasars and to evaluate the accuracy that can be realistically 
achieved on the value of Ho- 

Methods. We consider 14 lensed quasars that are candidates for time-delay monitoring and model them in detail using pixelized 
lens models. For each system, we provide a mass map, arrival-time surface and the distribution of predicted time-delays in 
a concordance cosmology, assuming Hq 1 = 14 Gyr (Ho = 70 in local units). Based on the predicted time-delays and on the 
observational circumstances, we rate each lens as 'excellent' or 'good' or 'unpromising' for time-delay monitoring. Finally, we 
analyze simulated time delays for the 11 lens rated excellent or good, and show that Ho can be recovered to a precision of 5%. 
Results. In combination with COSMOGRAIL paper I on the temporal sampling of lensed quasar light curves, the present work 
will help design monitoring campaigns of lensed quasars. 



Key words. Cosmology: distance scale, cosmological parametei 
1. Introduction 

Gravitational lensing of distant quasars is one of many 
possible routes to Hq. It has unique advantages. First, 
lensing depends on well-understood physics: gravitation. 
Second, time-delay observations require modest resources, 
hence low-demand telescopes can make a significant con- 
tribution. Third, the galaxy models used to convert 
time-delays into Hq have made considerable progress in 
the past decade. As a result, time-delay measurements 
have become an increasingly active research topic. 1 So 
far there are 12 secure time-delays (Table QJ, of which 
10 yield an estimate of the Hubble constant — the 
lens identification in PKS 1830-211 remains controversial 
QCourbin et al. 20021 IWinn et al. 2002jl . and the lensing 
galaxy HE 0435-122 may be anomalous (Kochanek 2005). 



1 According to ADS, the original paper by |Refsdal (1964)| 
pointing out the connection between gravitational lensing time- 
delays and Ho was cited on-average once every two years 
through the 1960s and 70s, whereas nowadays it is cited about 
once every two weeks. 



- Gravitational lensing - quasars: individual 

The dominant uncertainty in measuring Hq from lens- 
ing is the non-uniqueness of lens mass profiles that can 
reproduce the observables. Before the non-uniqueness of 
mass-models was widely appreciated, researchers would 
usually fit a single family of mass models to data, 
leading to over-optimistic error bars. Experimenting 
with different kinds of mass model for the same data 
pointed to much larger uncertainties (Schechter et al. 
1997, ISaha fe Williams 19971 Bernstein fe Fischer 19991 
Kect on et al. 200'Ojl . More recently, procedures involving 
sampling an ensemble of models according to some prior 
are being preferred, in order to derive a more use- 
ful picture of the uncertainties HWilliams fc Saha"2 000 
ISaha fc Williams 20041 |Oguri et al. 2004a| Jakobsson et 
al. 2005). A fair summary of current Ho results from lens- 
ing is that the error-bars are competitive on the young- 
Universe (i.e., high Hq) side but the old-Universe side 
needs improvement. 

Clearly to reach the 5-10% accuracy claimed 
by some other techniques (e.g., Frccdm an et al. 2 001 
|Spergel et al. 20031 ), more time-delays are needed. But to 
run monitoring campaigns efficiently, it is important to 
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have preliminary estimates for time-delays — witness the 
tenfold range in the known values in Table — as well 
as to identify the most promising systems to monitor. 
This paper supplies such information. For a sample of 14 
lenses we provide predicted time delays with uncertainties, 
a rating of prospects as 'excellent', 'good', or 'unpromis- 
ing' based on both models and the observational situation, 
and finally an estimate of the precision on Hq obtainable 
from these lenses. A companion paper by Eigenbrod et 
al. (2005; COSMOGRAIL I) is devoted to the determina- 
tion of the optimal strategy to use in order to measure 
time-delays (temporal sampling of the light curves, object 
visibility and variability, contamination by microlensing, 
etc). Together, these papers help design an observational 
campaign. 

An ideal time-delay lensing system has the following 
features: 1- bright optical images, 2- large angular image 
separations (>1"), 3- light path unperturbed by nearby 
structure, 4- known or easy-to- measure lens redshift zi cns . 
Our sample of 14 has been selected using these criteria 
as a guideline, though not a strict requirement. We have 
considered only objects for which the time-delay can be 
measured in the optical. 

The main results of this paper are in Sect. |3]and0J 
which present ensembles of models for the 14 individual 
systems and then estimate the precision to which H could 
be recovered from them. But before going into details of 
the models, it is useful to preview the results and compare 
them with measured systems. We do this in Sect. |2 

2. Comparing observed and predicted delays 

It is possible to make a preliminary prediction of time 
delays from image positions before any modelling, by re- 
calling the scales involved. 

In lensing theory, the geometric part of the time de- 
lay is of the order of the image-separation squared times 
DHq 1 , where D is the usual dimensionless distance fac- 
tor depending on cosmology. 2 The total time delay will be 



smaller but of the same order. Saha (2004) shows that the 
longest time delay can be expressed as 



(1) 



where where 9\ , 82 are the lens-centric distances (in ra- 
dians) of the first and last images to arrive 3 and ip is 
a dimensionless factor that ranges within about 0-2 for 
quadruples and 2-6 for doubles. The expression in square 
brackets in Eq. Q has the elegant interpretation of the 
fraction of the sky covered by the lens, times Hq 1 . 

We now define an 'astrometric time delay' Ai as t r0 m 
by taking Eq. and setting ip to a fiducial value of 1.5 



2 We refer all time-delay predictions in this paper to the con- 
cordance cosmology (S7 m = 0.3, JIa = 0.7) and H^ 1 = 14Gyr 
(or Hq = 70 in local units). 

3 In this section, in order to summarize the time-delays of 
many lenses, we will make the brutal simplification of neglect- 
ing the second and third images in quadruples. 



Table 1. The 12 time-delays measured so far, with la 
error bars. Lens redshifts in parenthesis are either photo- 
metric or based on absorption lines in the quasar images. 



Object 


Type 


•2-lens 


A^astrom 


Atobs 




B0218+357 


AD 


0.68 


10 


10 


± 




J0951+263 


ID 


(0.24) 


11 


16 


± 


2 C 


B1115+080 


IQ 


0.31 


24 


25 


± 


4 d > e 


B1600+434 


AD 


0.41 


35 


51 


± 


4 f 


B0435-122 


CQ 


0.46 


41 


14 


± 


I s 


B1830-211 


AD 


(0.89) 


42 


26 


4-5 
-4 


h 


B2149-274 


AD 


0.50 


59 


103 


± 


12' 


B1608+656 


IQ 


0.63 


60 


77 


± 


3 j 


B1520+530 


ID 


0.72 


92 


130 


± 


3 k 


J0911+055 


SQ 


0.77 


119 


146 


± 


8 1 


B1104-181 


AD 


0.73 


345 


161 


± 




B0957+561 


ID 


0.36 


536 


423 


± 


l n 



Cohen et al. 



2000)1 

Jakobsson et al. (2005) 



Barkana (1997)| 
^Kochanek et al. (2005)1 



Biggs et al. (1999)| 
Schechter et al. (1997)| 



IBurud et al. (2002a) | 
^Burud et al. (2002b)| 
^Ofek fc Maoz, D. (2003)] 



Burud et al. (2000)| 
Lovell et al. ( 1998)| 



IFassnacht et al. (2002)] 
]Hjorth et al. (2002) | ' 



Pscoz et al. (2001)| 




10 100 

astrometric delay (days) 



Fig. 1. Plot of Ai bs against Ai as t r om for the known time- 
delay systems, showing that known and prospective sys- 
tems can be easily compared. Squares denote quadruple 
systems, triangles are for doubles. 



for all quadruples and 4 for all doubles. This is a useful 
preliminary predictor of time delays, as we will see below. 

Table^gives the astrometric and actual observed time 
delays for the 11 known time-delay systems (disregarding 
here the middle two images in quadruples, i.e., images 2 
and 3 in the figures). The 'type' refers to the morpholog 
ical classification introduced in 
AD = axial double, ID 



Saha k Williams (2003) 



inclined double, CQ = 
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Table 2. Predicted time-delays (and la error bars) sorted 
by increasing astrometric delay, for objects with no mea- 
sured time-delay. Lens redshifts in parenthesis are either 
photometric or based on absorption lines in the quasar 
images. 



Object 


Type 


■2-lcns 


^^astrom 




B1422+231 


LQ 


0.34 


8 




J2026-453 


IQ 


(0.5) 


14 


15l 2 6 


J1155+634 


AD 


0.18 


19 


or + 8 

oo_ 10 


J0924+021 


IQ 


0.39 


19 


127J 


J1650+425 


ID 


(0.5) 


46 


54t? 3 


J1335+011 


AD 


0.44 


47 




J1355-225 


AD 


(0.70) 


68 




J1131-123 


LQ 


0.30 


69 


137 6 i39 


J2033-472 


IQ 


0.66 


70 


79+33 
' z -20 


B1030+074 


AD 


0.60 


75 


15312? 


B0909+532 


ID 


(0.83) 


90 


79+10 


B1009-025 


AD 


0.87 


98 




B0818+122 


ID 


0.39 


111 




J0903+502 


ID 


0.39 


122 






10 100 

astrometric delay (days) 

Fig. 2. Plot of A£p re d against Atastrom for the prospec- 
tive time-delay systems. Error bars are 68% confidence. 
Squares denote quadruple systems, triangles are for dou- 
bles. 



core quad, LQ = long-axis quad, SQ = short-axis quad, 
IQ = inclined quad. 

Frg-fflplots the data summarized in Tabled It is strik- 
ing that while Ai b s ranges over a factor of 40, it tracks 
A£ as trom to a factor of 2. 

Fig. |21 and Table |3 summarize our time-delay predic- 
tions. To make these predictions we used the PixeLens 
code l|Saha fc Williams 2004j) to generate an ensemble of 
200 models for each lens, leading to an ensemble of model 
time-delays, which we interpret as the probability distri- 
bution for the predicted time-delays. 




10 100 

observed delay (days) 



Fig. 3. Plot of Ai pre( j against At ^ s for the current time- 
delay systems. Again, squares are for quadruple systems, 
triangles for doubles. 



How reliable are the time-delay predictions? Pixellated 
models generically involve a choice of prior (also called 
secondary constraints); if the prior is too different from 
what lenses are really like then the results will be incorrect. 
Our prior is basically the PixeLens default; in detail, we 
assumed the following: 



1. 



In most cases we required the mass profile to be 
inversion-symmetric about the lens centre. But if the 
lensing galaxy appeared very asymmetric, or the im- 
age morphology was unusual, we let the mass profile 
be asymmetric. 

If there was evidence of external shear from the lens 
environment and/or the image morphology, we allowed 
the code to fit for constant external shear. That is to 
say, we allowed a contribution of the form 71 (8^ — 9y) + 
^li^x^y to the arrival time, with adjustable constants 
7i,72- 

The density gradient must point within 45° of the lens 
center (thus ensuring that the lens is centrally concen- 
trated). 

The radial mass profile must be steeper than 8~ 5 . 
That implies a 3D profile steeper than r -15 , which is 
consistent with available estimates from stellar or gas 
dynamics; for example, Binney et al. (1991) report an 



r -i.7S p ro fii e near the Galactic centre. 
5. The density on any pixel must be < twice the average 
of its neighbours, except for the central pixel, which 
can be arbitrarily dense. 

As a test we 'postdicted' the time-delays in the known 
systems. The results are summarized in Fig. We find 
that our prior tends to overestimate the time-delays for 
the systems with the largest angular separations, perhaps 
because these lenses have a significant cluster contribu- 
tion and the profiles are much shallower than in our prior. 
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One of the discrepant lenses is PKS 1830-211, which has 
a double lens galaxy. The two others are B0957+561 and 
J091 1+055, which both have significant contribution by 
a group or cluster of galaxies along the line of sight. But 
predicted time delays of less than 200 days appear reliable. 
The candidate lenses are all in the reliable regime. 

3. Individual systems 

We now proceed to discuss individual lenses, grouped by 
similar morphology. 

For each lens, we show three kinds of plot. First, there 
is a mass map of the ensemble-average model. The con- 
tours in the mass maps are in logarithmic steps, with each 
step corresponding to a factor of 10 4 (like a magnitude 
scale). The critical density contour is always the third 
from outermost. Second, we have plots showing saddle- 
point contours. These plots also show the source position 
in the ensemble-average model. The detailed placement of 
the saddle-point contours and the inferred source varies 
across the ensemble, but the qualitative features arc ro- 
bust. In particular, the saddle-point contours make the 
time-ordering of images obvious. We will refer to individ- 
ual images by their time order: 1,2 for doubles or 1,2,3,4 
for quadruples, meaning that the image labelled 1 varies 
first, then 2, etc. Third, we have histograms for the pre- 
dicted time delays between different pairs of images in 
each lens. 

After modelling each lens, we rate its prospects as a 
time-delay system as 'excellent', 'good', or 'unpromising', 
based on how well-constrained the time-delays are and on 
the comparative ease of monitoring and photometry. 

We remark that the modelling process really produces 
a predicted distribution for H At. In the present work 
we insert a fiducial value of Hq to obtain a distribution 
for At, but one can equally insert a measured value of 
At (if available) and obtain a distribution for Hq. But if 
two or more time delays become available for a quadruple, 
their ratio provides a new constraint on the lens, and the 
modelling code must be run again. 

3.1. Axial doubles 



J1155+634 [Fig.|U discovery: | Pindor et al. (2004)| The 

separation AO = 
galaxy is only ~ 
measurement z\ en 



1.83" is relatively large, but the lens 
0.2" from the fainter image. Also, the 
s = 0.1756 is somewhat insecure because 
the inferred galaxy absorption features are amongst the 
Lya forest lines. As a time-delay prospect, this system 
appears unpromising. 

J1355 225 [Fig-El discovery: |Morgan et al. (2003a)| also 
known as CTQ 327. The quasar images are bright and 
the angular separation is moderate: AO = 1.22". Models 
include external shear corresponding to further mass to 
the NW or SE. We rate this system as a good time-delay 
prospect. 

J1335+011 [Fig. E| discovery: |Oguri et al. (2004b)| We 




20 40 
Predicted time delay (days) 



Fig. 4. Models of J1155+634 (axial double). See text in 
section El for the format. Prospects: unpromising. 





40 60 80 100 120 
Predicted time delay (days) 



Fig. 5. 

good. 



Models of J1355-225 (axial double). Prospects: 




20 40 60 

Predicted time delay (days) 



rate this system as an excellent time-delay prospect. 



Fig. 6. Models of J1335+011 (axial double), 
excellent. 



Prospects: 
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Predicted time delay (days) 



20 40 60 80 

Predicted time delay (days) 



Fig. 7. Models of B1030+074 (axial double). Prospects: Fig. 9. Models of J1650+425 (inclined double). Prospects: 
unpromising. excellent. 






100 200 
Predicted time delay (days) 






20 40 60 80 100 

Predicted time delay (days) 



Fig. 8. Models of B1009-025 (axial double). Prospects: 
good. 



Fig. 10. Models of B0909+532 (inclined double). 
Prospects: excellent 



3.2. Inclined doubles 



B1030+074 [Fig. |7] discovery: Xanthopoulos et al. 
(1998). Like J1155+634 it has a relatively wide separa- 
tion but a second image is faint and very close to the 
galaxy. There is evidence for variability. The peak in the 
predicted time delays near 180 days is interesting, but it is 
probably not wise to over-interpret, given the resolution 
of the models used in this paper. Because of the difficulty 
of accurate photometry on the second image, we rate this 
system as unpromising. 

B1009 025 [Fig. discovery: |Surdej et al. (1993)| Its 



clean morphology, evidence of variability and a nearby 
foreground QSO usable as a standard PSF all make this an 
attractive target. However, the combination of an approxi- 
mately half-year time-delay and a near-equatorial location 
is awkward (see Eigenbrod et al. 2005 for more details). 
We rate time-delay prospects as good. 



J1650+425 [Fig. OH discovery: |Morgan et al. (2003b)| It 

has significant external shear, probably from a group 
galaxy to the E. The high declination of the objects makes 
it possible to observe it almost continuously from the 
northern hemisphere. This system is an excellent time- 
delay prospect. 

B0909+532 [Fig. EH] discovery as multiply imaged: 
Koc hanek et al. (1997)| The lensing galaxy is very faint, 



which caused some early controversy until the issue was 



settled by Oscoz et al. (1997) and Lubin et al. (2000) 



The morphology and models indicate significant external 
shear from mass to the NE or SW, but the galaxies respon- 
sible have not yet been identified. Both quasar images are 
very bright, and their separation is A9 = 1.17". Its zi en s 
is still insecure, but assuming that problem is solved, this 
system is an excellent prospect. 
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1 


1-2 


nr 






■ , , r, , 







20 40 60 80 100 120 140 160 
Predicted time delay (days) 



Fig. 11. Models of B0818+122 (inclined double). 
Prospects: good 




20 40 60 80 100 120 140 160 180 
Predicted time delay (days) 



Fig. 12. Models of J0903+502 (inclined double). 
Prospects: good. 



B0818+122 [Fig.ITI] discovery: Hagen & Reimers (2000). 
A chain of galaxies to the NE contribute a large external 
shear. The fainter image is very close to the main lensing 
galaxy, and about the same brightness. Overall, prospects 
appear good. 



J0903+502 [Fig. U2 discovery: |Johnston et al. (2003) 
There are several group galaxies in addition to the main 
lensing galaxy, with one galaxy to the SW probably the 
major contributor of external shear. Both quasar images 
are faint, R ~ 19 — 20, so monitoring is difficult with a 1 m- 
class telescope. The distribution of predicted time-delays 
is narrow. Overall, we rate prospects as good. 

3.3. Long-axis quadruples 



B1422+231 [Fig. EI discovery: |Patnaik et al. (1992)| It 

is a radio emitter with extremely accurate image posi- 
tions. There is evidence of variability. Strong external 




1.0 2.0 3.0 4.0 

Predicted time delay (days) 




10 20 
Predicted time delay (days) 



Fig. 13. Models of B1422+231 (long-axis quadruple). 
Prospects: good. 



shear comes from a galaxy group to the SE. Time-delays 
between the close triplet of images may be too short to 
measure in the optical, but the delay to the fourth image 
can be expected to be useful. The lensing galaxy is compa- 
rable in brightness to the faint fourth image, which com- 
plicates the photometry. Overall, prospects appear good. 

It is 



J1131 123 [Fig. ED discovery: |Sluse et al. (2003) 
a quadruple with large separation: 3.69". It is very like a 
larger sibling of B1422+231. Morphology and models indi- 
cate significant external shear from mass to the WNW or 
ESE. There is evidence for intrinsic variability. Structures 
in the Einstein ring arc likely to offer additional model 
constraints, though they also contaminate the flux from 
the faint fourth image. Overall, prospects appear excel- 
lent. 

3.4. Inclined quadruples 

J2026-453 [Fig.lTSl and J2033-472 [Fig. Q3J discovery: 
Morgan et al. (2004) In J2026-453, morphology and mod- 
els indicate external shear from mass to the E or W. This 
system so far lacks a zi cns ; we assumed 0.5, which is plau- 
sible given the colours of the galaxy. The morphology of 
J2033-472 suggests an asymmetric lens, and accordingly 
we have considered asymmetric models. We rate time- 
delays prospects as good for J2026 and excellent for J2033. 



J0924+021 [Fig. 02] discovery: |lnada et al. (2003)] This 
is a complex and evidently asymmetric lens, but can be 
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2 4 6 

Predicted time delay (days) 




100 200 300 

Predicted time delay (days) 



Fig. 14. Models of J1131-123 (long-axis quadruple). 
Prospects: excellent. 






6 8 10 12 
Predicted time delay (days) 




4 6 
Predicted time delay (days) 



Fig. 15. Models 
Prospects: good. 



of J2026-453 (inclined quadruple). 




20 40 60 80 100 

Predicted time delay (days) 




2 4 6 8 10 

Predicted time delay (days) 




20 40 60 80 100 

Predicted time delay (days) 

Fig. 16. Models of J2033-472 (inclined quadruple). 
Prospects: excellent. 

well-modelled and leads to relatively tight predictions for 
two of the time delays. However, image 3 is very faint, 
which Keeton et al. (2005) argue is the result of microlens- 
ing. This greatly complicates the measurements of time 
delays, so we current rate this lens as an unpromising 
time-delay prospect. 

4. Predicted precision for the Hubble time 

In the previous section, after considering detailed mod- 
els as well as observational circumstances of all 15 
lenses, we concluded that 5 systems are excellent candi- 
dates for time-delay monitoring (J1650+425, J2033-472, 
B0909+532, J1335+011, J1131-123), and 6 are good can- 
didates (J1355-225, J0903+502, B0818+122, B1009-025, 
B1422+231, J2026-453). We now ask how accurately H^ 1 
can be inferred if the 5 excellent candidates, or if the 11 
excellent or good candidates, have their time delays mea- 
sured accurate to 1 d. We do not expect that either see- 
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Hubble constant (local units) 
80 70 60 50 40 



10 20 
Predicted time delay (days) 




2 4 6 8 10 12 14 
Predicted time delay (days) 



Fig. 17. Models of J0924+021 (inclined quadruple). 
Prospects: currently unpromising. 



nario will be what transpires in the future. We expect that 
some of these 1 1 lenses will yield accurate time delays over 
the next 2-3 years, while some existing time-delay mea- 
surements are refined. But the 5-lens and 11-lens cases are 
reasonable surrogates for a future set of available measure- 
ments. 

Fig. 1181 shows the recovered Hq 1 from simulated time 
delays of the 5 excellent candidates. For each lens we 
took a random model (from the ensemble of 200), read 
off its time delays rounded to the nearest day, and then 
took them as simulated time delays. Any model delays 
of < Id we treated as unmeasured. Using PixeLens, we 
then modelled the 5 lenses simultaneously from the ac- 
tual image positions and these simulated time delays. The 
model-ensemble had 200 members, each member consist- 
ing of models for all 5 lenses sharing a common -ff^ 1 
IjSaha fc Williams 2004|l . Fig. HU shows the resulting 200 
values of Hq 1 after binning. We see that the original in- 
put value 14 Gyr is recovered with < 10% uncertainty at 
68% confidence, and no discernable bias. Also, the uncer- 
tainties are asymmetric. 

Fig. 1191 shows the result of a similar exercise using the 
6 good candidates. The uncertainties are somewhat larger 
than in Fig. 1181 and similarly asymmetric. Also, there is a 
bias, in the sense that the median value is 15.0 Gyr rather 
than 14.0 Gyr; but the bias is in the 68% confidence range 
and hence not significant. 




2(1 

Hubble time (Gyr) 



Fig. 18. Hubble constant/time recovered from simulated 
time delays of the 5 excellent candidates. The confidence 
intervals (read off by sorting the unbinned values) are 
M.Otj^Gyr at 68% and 14.0t^Gyr at 90%. 



Hubble constant (local units) 
80 70 60 50 




20 

Hubble time (Gyr) 



Fig. 19. Hubble const ant /time from simulated time de- 
lays of the 6 good candidates. The confidence intervals 
are 15.0±^ Gyr at 68% and lS.Ol^ Gyr at 90%. 



Hubble constant (local units) 
80 70 60 




14 16 
Hubble time (Gyr) 

Fig. 20. Hubble constant /time from simulated time de- 
lays of the 11 excellent or good lenses. This histogram 
is simply the product of the two previous ones. From 
a narrower binning, we read off the confidence intervals 
14.2±°;£Gyr at 68% and 14.4±J jGyr at 90%. 

The 5-lens and 6-lens ensembles just described are 
independent. Hence we can simply multiply their his- 
tograms. Fig. 1201 shows the result. We recover Hq 1 with 
an uncertainty of about 5%, at 68% confidence. 

These results show that the Hubble time can be recov- 
ered to 5% precision even if we allow for a large diversity 
in possible mass distributions (or prior). But there is a 
caveat, which needs to be addressed before a claim of 5% 
accuracy (rather than precision) can be made. Currently, 
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a mass distribution that satisfies the lensing constraints is 
either allowed by the prior as a plausible galaxy lens, or 
rejected; there is no weighting in the prior. Properly, the 
prior should weight mass models according to their abun- 
dance in the real world of galaxies. Lack of weighting will 
introduce a bias. (This prior-induced bias is different from 
the small statistical bias seen in Fig. EH) The blind tests 
in Williams & Saha (2000) would have detected biases if 
they were around 20% or more. But prior-induced biases 
at the 5% level remain untested for. Finding them and 
then eliminating them through a weighted prior could be 
done by calibrating against galaxy- formation models, and 
is an essential theoretical program needed to complement 
the observations. 



5. Discussion 

In this paper we do three things: first, we introduce a 
simple rough predictor for time delays At as t rom j second, 
we make model predictions for 23 time delays covering 14 
lenses, and finally we estimate the precision in the Hubble 
time inferred from simulated data on the 11 best lenses. 
The main conclusion is that no single lens can usefully 
constrain Ho, but time delays accurate to ~ Id on > 10 
lenses can yield Hq accurate to 5%. 

In the histograms in Figs. ITHTTl typically 90% of the 
area ranges over a factor of two in time delays. Hence, a 
monitoring program can have 90% confidence in succeed- 
ing — provided the quasar is sufficiently variable — if the 
sampling allows for the appropriate 90%-range of possible 
time delays. There is no single characteristic shape for the 
histograms, but the pattern of a low-end tail and a high- 
end cliff is common. The large uncertainty in the predicted 
time delays reflects the large variety of mass models mod- 
els that can reproduce the observed image positions in any 
given lens. The prior we have for deciding what is an allow- 
able mass model for a galaxy is very conservative, so our 
models have more variety than reality. But not very much 
more — if all real galaxy lenses belonged to some known 
parametrization, then error-bars on Hq from fitting such 
parametric models to observed time-delay systems would 



overlap, but as is evident from Fig. 12 in Courbin (2003) 
those error-bars do not overlap. 

Actually, the basic results about predicted time delays 
are already present in the summaries Figs. ^ and [21 T° 
interpret these figures, recall that Ai as t r om makes a pre- 
liminary prediction for the time-delay, while the deviation 
from the oblique line depends on the details of mass dis- 
tribution and lens morphology. From Fig. we see that 
Aiastrom gets us to within a factor of two of the observed 
values. Now, the error bars — note that the error bars in 
tables and figures are 68% confidence — in Fig. [5] are gen- 
erally less than a factor of two; thus detailed modelling 
does provide a better prediction than Ai as t r0 m alone, but 
not dramatically better. We can also see that several of 
the error bars in Fig. |3 are shorter on top, thus indicating 
a low-end tail and a high-end cliff. 



The image morphology of a lens is correlated with the 
uncertainty in the time delays, especially in quadruples. 
Core quadruples generically have short time delays and are 
unlikely to be useful for time delays; the case of B0435- 
122 is illustrative. Long- and short-axis quadruples, hav- 
ing three images close together, are likely to have only 
one measurable time delay. Inclined quadruples are the 
most promising, since they usually have two time delays 
in the measurable range, and sometimes three. Among 
doubles, inclined systems tend to be somewhat better 
constrained than axial systems. Significant asymmetry in 
the lens is a disadvantage, but in compensation, asym- 
metry increases the chance of having three measurable 
time delays. Thus B1608+656 is an asymmetric inclined 
quadruple with three measured time delays; J2033-472 
may prove to be another, and is among our excellent can- 
didates. Surprisingly, a large external shear appears to 
reduce uncertainties in the time delay. This is particularly 
noticeable in the inclineddoubles J1650+425, B0909+532, 
and J0903+502, and the short-axis quads B1422+231 and 
J1131-123. The reason is not clear; it may be that since 
external shear reduces amount of mass needed in the main 
lens to produce multiple images, it reduces the available 
model-space. 

Finally, the results from combining several lenses are 
very encouraging. Assuming time delays accurate to 1 d we 
find that the model-dependent uncertainty in Hq 1 reduces 
to less than 10% on combining the 5 best lenses, and about 
5% on combining the best 11 lenses. The uncertainties are 
asymmetric, with the lower limit on the Hubble time being 
tighter than the upper limit. More work needs to be done 
on the model prior before we can truly attain 5% accuracy, 
but meanwhile our results help provide both motivation 
and observing strategies for accurate time-delay measure- 
ments. 
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